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Abstract 

The recent data for E866 on the xp dependence for charmonium suppression 

in pA colhsions at 800 GeV are analyzed using a time- and energy-dependent 

preformed charmonium absorption cross section afy^^ (r, yjs^^) . For ^/s = 

10 GeV the initiahy (r = 0) produced premeson has an absorption cross 

section of Cpr — 3 mb. At the same energy but for r ^ oo one deduces for 

the total cross sections (rfl^^ = (2.8 it 0.3) mb, a%l^ = (10.5 it 3.6) mb. The 

date are compatible with a formation time = 0-6 fm/c. 
PACS: 24.85. =p, 12.38.Mh, 25.75.Dw 
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Typeset using REVT^ 



The formation of a charmonium meson J / ip or ip' from the originally produced cc pair via 
gluon fusion is a point of discussion and disagreement in the community which investigates 
the charmonium suppression in nucleus-nucleus collisions What is the nature (color, 

size) of the premeson when it travels through nuclear matter? How long does it take to be 
a fully developed J/ip or ip'^ Are the formation times different for J/ip and ip', the ip' being 
considerably larger than the J/V"? What are the total cross sections (jfot'^ ■, '^tot^ j for ths 
fully developed charmonia on a nucleon? These questions can be possibly answered using 
the recent and still preliminary data of the E866 collaboration 0. The analysis will be 
given in this paper. Since the data are still preliminary, this letter stresses the method of 
the analysis, while the numerical results may still undergo certain modifications. 

In an experiment pA — > ipX at 800 GeV, where ip stands for J/ip and the suppression 

^^^^^^ = Ada{pN J/^X) 

has been measured as a function oi xp for the produced ip in the range —0.1. We will 
concentrate in this paper on the interval —0.1 xp 0.25 where formation time effects are 
expected to be of particular importance and where effects of the production time (coherence 
time) can be neglected. The effects of the formation time are determined by the Lorentz 
factor 7(xir) of the charmonium with respect to the target nucleus. In the experiment pA 
at 800 GeV, 7 varies between a value of 7 a,t xp = —0.15 and a value of 47 at xp = +0.15. 
For a hypothetical value of the charmonium formation time Tf = 0.5 fm/c in the c.m.s. and 
for a premeson which is produced in the middle of a nucleus with radius Ra = 6 fm, a ip 
observed with xp = —0.15 is essentially fully developed inside the nucleus and absorbed 
with the asymptotic total cross section af^f , while ior xp = 0.15 the ip traverses the nucleus 
in the form of a premeson. 

On the other hand, the coherence time tc of charmonium production, which is the 
quantum-mechanical uncertainty in the time of production of the premeson is quite short 
in the interval of consideration. In the light-cone approach tc can be treated as the lifetime 
of a fluctuation containing the cc. It is given by the energy denominator corresponding to 



such a fluctuation, 

t - ^ (2) 

where is the effective mass of the fluctuation. We assume here that the cc pair is 
produced via gluon-gluon fusion accompanied by gluon radiation which is needed in order 
to end up with a color neutral cc pair. It would be incorrect to say that the gluon is radiated 
after its interaction with the target. As usual with bremsstrahlung it takes place both before 
and after the interaction. This is a quantum-mechanical uncertainty. In the light-cone 
approach tc corresponds to the lifetime of the whole fluctuation G ccG. This is why 
Eq. (0) contains the gluon energy Eg in the numerator and in the denominator the effective 
mass of the qqC , 

= — ^ + ^ . (3) 

a{l — a) 1 — a 

Here kr and a are the transverse momentum and the fraction of the initial light-cone mo- 
mentum carried by the radiated gluon, respectively. Mgc is the effective mass of the cc pair 
which we assume to be of the order of the charmonium mass. It is important to note that 

cannot be small since the radiated gluon has to resolve the inner structure of the cc pair 
in order to make it colorless. Therefore, it should be larger than the inverse cc separation, 
kT > 1/r^ ~ rric- At the same time, the fraction a should not be large, otherwise it will 
produce a shift in the initial gluon momentum, which falls off steeply like (1 —xi)^, resulting 
in a strong suppression. A detailed calculation leads to values tc = 0.25 fm/c at xp = —0.15 
and tc = 4.1 fm/c at xp = 0.15 in the lab system and Tc = 0.04 fm/c and Tc = 0.08 fm/c in 
the premeson c.m.s., respectively. Indeed we will see ^ rj. For values of 0.3 < xp < 0.8 
one has tc > R and the formation time concept becomes doubtful, another mechanism seems 
to be responsible for the observed suppression. 

In the Xp interval [—0.1, 0.25] we model the suppression by assuming that the premeson is 
created instantaneously and is absorbed according to a time and energy dependent effective 
absorption cross section cr^^{T, ^/s^n) for an experiment, in which a charmonium ip is 
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observed. (The effects of the decays J/^J and ip' ^ J/ip will be discussed later.) The 

classical expression 



S^{xf) = -^j d% J dz pA{b, z) exp I- dz' pA{b, z')a^^^ (7(3?^' V^'^^^^^))} ' *^4) 

—00 

is used with the nuclear density distribution p^ih^z) normalized to A. This expression is 
an approximation to exact solutions for the evolution of a wave packet with charmonium 
quantum numbers propagating through a medium treated in the quark [Q or in hadronic 
||T0| representations (see also below). We believe that Eq. (|D is more intuitive and easier 
to use than the exact ones. Our aim is to determine the dependence of cr^^f (r, ^s^n) and 
'^abs(''"' v^s^at) from a fit of Eq. @) to the E866 data. 

We use the following expression for o-^j^(r, y/s^) with two adjustable parameters Eg 
and Soo which correspond to the effective absorptive cross section at short and long times, 
respectively, 

^tZ (r, v^) = Poo + (Soo - So) cos(AM r)] (f^) • (5) 
The dependence on the energy y/s^N is deduced from photoproduction experiments {•yp — >■ 



J/ipp) with A = 0.4 |n| The form of the dependence on the time r is derived within the 
following quantum mechanical model of two coupled channels for the evolution of a color 
neutral cc pair [QQ: The time dependent premeson state \cc{t)) with the quantum numbers 
of J/ ip and ip' can be expanded in a complete set of hadronic states of which we keep only 
the lowest two ones, J/ip and and we may use spinor representation: We denote by 



|cc(0)) 



(6) 



the initially produced superposition of the J/ip and The |cc) wave packet in an interacting 
environment is described by the equation. 



^ We assume that the state |cc) includes all the Fock components with additional gluons and sea 
quarks, |cc) = \cc)o + \ccG) + ... 



\dcc{T)) 
dr 



(7) 



with 
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(8) 



O"00 O"oi 

where Q is the mass matrix and T is the interaction amphtude operator containing diagonal 



and off-diagonal amplitudes, Uoo = {J / '>p\d'\J / Tp) , ctqi = {J/M^W)^ ^^c. |]T0 



For a |cc) created at the point {h,z) in the nucleus and observed asymptotically as a -i/^ 
one has the transition probability 



cc 



'■-)) 



(9) 



mc-cmv 

where |cc(^^)) (it depends also on h) is the solution of Eq. (|^ with the initial state (|) at 
the point with coordinates {h,z). 

Expanding expression (^ in up to the first order we get, 



with AM = M^/ - Mji^. 



dz'pA^b, z) 
dz'pA{h, z) 



(Too + -RcTio COS AM 



z — z 
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1 f,-^^---- 



0-11 + -^c^oi COS I^AM 
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(10) 



The expressions in square brackets in Eq. (0) are the time dependent effective absorption 
cross sections and are of the form assumed in Eq. (^. These effective cross sections may be 
positive and negative. In the latter case one observes an enhanced production if one uses a 
nuclear target 

We have calculated the suppression function S^{xf) in Eq. (H) for Tungsten A = 182 
and Beryllium A = 9 with a uniform density, and R = tqA^^^ with tq = 1.14 fm. The two 
parameters Eq, and Sqo for each species of charmonium have been determined by fitting 
Eq. ([H) to the data of J/ip and ip' suppression, respectively. We have used MINUIT-Hesse 
from CERNLIB. Fig. [l| shows the fits. The numerical values of the parameters together with 
their errors and the values Xdof given by the fit routine are displayed in Table 1. 



In the following discussion of the results we always take ^/s^ = 10 GeV. 

For T — > oo, the oscillating term in the parameterizations Eq. (j^) does not contribute 
in the integral for the suppression. This is the situation of the fully developed charmonium 
and the parameter Ego can be identified with the total cross sections, 

al^J^""" = (5.0 ± 0.4) mb, atl" = (10.5 ± 3.6) mb. (11) 

We have set J / ip in quotation marks, since the total cross section a^oJ^ ^ is the one for a 
situation, where the observed J/ip originates with probability pi ~ 0.6 from directly formed 
J/ip, and probability p2 — 0.3 and ps ~ 0.1 from the decay of Xc and ip', respectively. Thus 
the total effective absorption cross section seen in the ^^J/ip" channel is a superposition of 
the contributions of J/ip, Xc, and tp'. If we correct for this effect, assuming that total cross 
sections are proportional to (r^) one has 



J/tpN _ "J/4!"N 
^tot — ^tot 



Pi + 



= 2.8 ± 0.3 mb , (12) 

where we have used {r'^fj^^ = 0.42 fm, (r^)!/^ = 0.67 fm, {r%i'^ = 0.85 fm jg. The value 
in Eq. (|12|) can be compared with afl^^^ = (3.5 ± 0.7) mb obtained from an analysis of 
photoproduction data 7p J/ipp using the modified vector dominance model [^. 
The ratio of the values from Eqs. (pUD and (|T^) gives 

-^ = 3.8 ±1.3, (13) 
which is close to the value derived from ip' photoproduction |jll| and coincides with the 



prediction based on the ratio of the corresponding values of (r^) which leads to 4 when the 



calculations by Buchmiiller |]T2[ are used. 

Next we discuss the size of the absorption cross sections for a premeson directly after 

creation, i.e. at r = 0. Then according to Eq. (|^), the cross sections are related to Sq. 
Using the data from Table 1 we have 

2.7 ± 0.1 mb J /%p observed 

3.8 ± 0.6 mb ip' observed . 
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The two values are rather small compared to the asymptotic cross sections, supporting the 
idea that the initial cc system is rather small. The values in Eq. (|T^ are equal within the 
large error bars, although there is no compelling reason why they should be the same. For 
instance, one of the origins for a difference is the admixture of the Xc-states in the J/ij) 
channel, which at r ^ cannot be corrected for easily since the (r^) law is not applicable. 

Within the two channel model the time dependence of the effective absorption cross 
section, Eq. (|^) is given by the characteristic time 1/AM = 0.3 fm/c, where AM is the 
mass difference between the states J/ip and ip' . We have not varied this characteristic time 
since only few data points are available with rather larger error bars. However, since the 
parameters Eq and Sqo deduced from the experiment yield a consistent picture and values of 
X^dof < -'-) time dependence derived from the two channel model may not be unreasonable. 
The oscillating dependence of the cross section in Eq. makes the physical interpretation 
of the characteristic time somewhat difficult. Since the cross sections always appear under 
the integral we may define a formation time ri/2 by the requirement that the contribution 
of the oscillating term is reduced by 50%, i.e. we define ri/2 by 

^ ^ drcos AMr = i , (14) 



Ti/2 Jo 2 

which happens at ri/2 = 1.9/AM = 0.6 fm/c. 

The previous analysis is based exclusively on the data of E866 and has yielded the time 
and energy dependent absorption cross sections o"^bI(T, ^/s^n) parametrized as in Eq. @. 
We have checked, whether our results are consistent with other data, for instance, measure- 
ments of J/ip suppression at 200 GeV (NA3) and ratios of ip' to J/ip suppression at 450 
and 200 GeV (NA38/50). In the experiment of NA3 [|l|], where Pt and have been the 



targets, only four data points are available for xp ^ 0.4, where formation time effects are 
supposed to be the dominant mechanism. We have varied only Ego, while setting Sq at the 
value derived from the 800 GeV data. The results from the fit are also given in the Table 1 
and agree with the results from E866. 

Furthermore, we have used the absorption cross sections as deduced from E866 to calcu- 



late the ratio of ip' /J/ip production on nuclear targets which has been measured by NA38 
and NA50 for several nuclei and energies of 200 and 450 GeV. Results are shown in Fig. 2. 
Within the uncertainties induced by error bars of the fitted parameters, the data on / J /if) 
suppression can be understood with the time dependent absorption cross section. 

We summarize: The data of the E866 experiment for the xp dependence of the J/%Ij 
and ip' suppressions (—0.1 < Xi? < 0.25) have been analyzed with the hypothesis that the 
variation is entirely due to a time- and energy- dependent effective absorption cross section 
for a premeson which evolves in time during its passage through the nucleus. The analysis is 
based on a set of hypotheses and approximations which we list, but whose accuracies cannot 
be estimated: 

• We assume that formation time effects dominate the XiT'-dependence of nuclear sup- 
pression for xp ^ 0.25 where the coherence time is rather short. 

• We assume that the formation times for J /ip and if)' are equal. It is frequently assumed 
based on the picture of classical expansion that the ip' needs longer time to form than 
the J/ip because of its larger radius. This might be not true in quantum mechanics. 
For instance, for the oscillator potential the formation time is size independent (the 
period of oscillation of a pendulum is independent of the amplitude). 

• The energy and time dependence of the absorption cross section factorizes 
o"abs(T, -\/s) = o"(t)(s/so)'^. The dependence of A on the transverse cc separation is 



rather weak ||15[ and we fix it at A = 0.2 which follows from J/\l/ photoproduction 
data. 

We have assumed one exponential in each channel when calculating the suppression 
function S^{xf), although the situation is more complicated (coupled system of J/ip 
and ip' , contribution of Xc eic.) as is suggested by the evolution equation (^ for two 
coupled channels and the exact solution (P) which can be found in |]10 . 



Despite these uncertainties, the analysis has provided a coherent description of the data from 



E866, NA3 and NA38/50. We have extracted values for a^l^^ and 0"^/^, which agree with 
values extracted from photoproduction experiments and follow the systematics of values of 
hN total cross sections a'^^ as a function of {r'^)h- For small times, r — >■ 0, one finds smaller 
absorption cross sections, as expected, if the premeson is small in size. No statement can be 
made about the color structure. Although we have not varied the formation time, the fit is 
very good for the formation time ri/2 — 0.6 fm/c derived from the two channel model. 

After our analysis was complete Ref. |TB| has appeared in which the same data are 
analysed as in this work but with a quite different approach. It covers the whole range of 
Xp, and treats the evolution of the cc pair purely classically. 
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FIG. 1. Data from E866 for the ratio of charmonium production in pW and pBe collisions 
together with our best fits using Eq. (4). 
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FIG. 2. Results for the ratio B'^i^a{ip')/Bi^na{J/ip) calculated from fitted parameters in Table 1 
are compared with pA data at and 200 and 450 GeV (data from NA51 and NA38). The solid curve 
corresponds to the calculation with parameters at their central values, while the upper and lower 
curves characterize the uncertainties induced by the error bars of the fitted parameters. 
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TABLES 





(exp.) 


Soo[mb] 


So[mb] 


Xdof 




(E866) 


5.0 ±0.4 


2.7 ±0.1 


0.84 




(E866) 


10.5 ±3.6 


3.8 ±0.6 


0.28 


J/i; 


(NA3) 


6.8 ± 1.7 


[2.7] 


1.6 



TABLE L Values for the parameters Ti^o and T,q in the parametrization of the absorption cross 
section Eq. @ as obtained from the least square fit to the data. For the NA3 experiment only the 
parameter Sqo has been fitted, while the values of Sq has been set at 2.7. 
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